Cosmic Ray Energy Spectrum from Measurements of Air Showers 
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This review focuses on high-energy cosmic rays in the PeV energy range and above. Of particular 
interest is the knee of the spectrum around 3 PeV and the transition from cosmic rays of Galactic 
origin to particles from extra-galactic sources. Our goal is to establish a baseline spectrum from 10 14 
to 10 20 eV by combining the results of many measurements at different energies. In combination with 
measurements of the nuclear composition of the primaries, the shape of the energy spectrum places 
constraints on the number and spectra of sources that may contribute to the observed spectrum. 
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I. INTRODUCTION 

In the 100 years since Victor Hess's balloon night that 
marks the discovery of cosmic rays, a great deal has 
been learned about their composition, propagation and 
sources. We know that most particles originate from 
sources in the local galaxy, having spent on average 10 7 
years in diffusive motion in the interstellar medium (ISM) 
and the galactic halo before being lost to intergalactic 
space [l| . We know that the power requirement to main- 
tain the intensity of cosmic rays at a near constant value 
is a few per cent of the energy available in supernova 
explosions and that there is a well developed theory of 
non-linear, diffusive acceleration by shocks driven by ex- 
panding supernova remnants (SNR) [2J. We know that 
diffusion in the interstellar medium depends on energy in 
such a way that higher energy nuclei escape more quickly 
than lower energy ones, but the shape of this energy de- 
pendence is still uncertain at high energy Q. 

Measurements of diffuse gamma radiation from < 1 
to > 100 GeV by the Fermi Satellite [4j confirm this 
general picture of the origin and propagation of galactic 
cosmic rays. In the TeV range, ground based gamma-ray 
telescopes have observed emission of >TeV gamma-rays 
from supernova remnants, as expected if there is parti- 
cle acceleration to energies of > 10 TeV However, 
because electrons are more efficient radiators than pro- 
tons, the relation of these observations to the fluxes of 
high-energy protons and nuclei observed at Earth is not 
completely clear. In particular, the maximum energy to 
which protons and nuclei can be accelerated by supernova 
remnants is not yet directly established. 

The picture outlined above is driven by direct mea- 
surements of cosmic rays with detectors near the top 
or outside of the atmosphere. Measurements with spec- 
trometers in high-altitude balloons and satellites mea- 
sure nuclear and isotopic composition in the GeV en- 
ergy range and spectra of individual elements up to 
TeV/nucleon [6|. Instruments with calorimeters and 
transition radiation detectors extend the energy measure- 
ments to ~ 100 TeV per nucleus. Recent measurements 
by CREAM and PAMELA @ support earlier obser- 
vations by JACEE [l(J [Tl| and ATIC [lj] that the spec- 



trum of helium is somewhat harder than that of protons. 
The PAMELA measurements show that the spectrum of 
both protons and helium become harder above a rigidity 
of 200 GV. The CREAM data suggest that this is also 
the case for heavier nuclei Q. Further new data are ex- 
pected soon from the AMS-2 detector that by now has 
worked on the International Space Station for more than 
one year. 

The cosmic-ray flux above 100 TeV amounts to about 5 
particles per square meter per steradian per day. There- 
fore, direct measurement of the cosmic-ray spectrum 
above this energy with balloon borne or satellite detec- 
tors is difficult. The highest direct measurement so far is 
still the series of measurements by Grigorov et at carried 
on Proton satellites [l3|, which extends to ~PeV. Data 
for the spectrum above a few hundred TeV comes from 
large air shower arrays on the ground. The air-shower 
experiments observe the cascades of secondary particles 
in the atmosphere initiated by the interaction of the high 
energy primary particles. In such indirect measurements, 
the information about composition is limited (at best) to 
determining the relative abundances of the main groups 
of nuclei. Moreover, the relation between the observed 
signal at the ground and the primary energy depends on 
the mass of the primary nucleus. In addition, the rela- 
tion between the observed signal and the primary energy 
depends on the model(s) of hadronic interactions used 
to interpret the cascades. Given the large uncertain- 
ties in energy assignment, it is not surprising that dif- 
ferent measurements give different results in the same or 
overlapping energy regions. Another complexity is that 
techniques for measuring the shower at the ground differ 
from one experiment to the next (for example, scintilla- 
tors compared to water Cherenkov tanks) . In this review 
we therefore use observed features in the energy spectra 
to cross-calibrate results from different arrays. 

Since air shower measurements are calorimetric in na- 
ture, the natural energy variable to use is total energy per 
nucleus. This contrasts with the direct measurements, 
which generally use energy per nucleon. The latter is 
the natural variable for studying propagation of nuclei 
because energy per nucleon is conserved in spallation. A 
variable different from both of these is relevant for in- 
terpreting observed spectra. This is magnetic rigidity, 
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FIG. 1: All particle cosmic ray spectrum from air shower experiments. (References in text.) 



defined as 

Pc 

R > £. (1) 

where P is the total momentum of a nucleus and Ze its 
electrical charge. Particles with the same rigidity and 
injection vector follow identical trajectories in a given 
magnetic field configuration. Rigidity is therefore the 
appropriate variable for interpreting changes in spectrum 
due to propagation and acceleration in magnetic fields. 
In particular, as first pointed out by Peters if there is 
a maximum energy to which protons can be accelerated 
in a source, then the protons will cutoff first, followed by 
helium, carbon, . . . according to 

E max {Z) = ZexR c = Zx E max {Z = 1). (2) 

The outline of the paper is as follows. We first describe 
briefly the different types of air shower experiments and 
summarize the data from each of the selected measure- 
ments. We then use features in the energy spectrum ob- 
served by different experiments to construct a tentative 
all-particle energy spectrum from 10 14 to 10 20 eV. The 
most prominent features are the knee around 3 x 10 eV 
and the ankle around 10 19 eV, both prominent in Fig. [T] 
As justification for this approach, we note that, to a high 
degree of accuracy, there must be a single spectrum at 



Earth. In the third major section of the paper we use 
this constructed all-particle spectrum as a template for 
discussing measurements of composition and possible im- 
plications for different sources. In this section we will as- 
sume the validity of the Peters cycle as written in Eq. [5] 
as a constraint on the energy dependence of different nu- 
clear components. We describe two fits to the data, each 
of which has three populations of particles with contrast- 
ing assumptions about the rigidity cutoff for each popu- 
lation. The first two populations represent cosmic rays 
from galactic sources and the third population is an ex- 
tragalactic component. Each population contains sev- 
eral groups of nuclei with assumed spectral indices as 
adjustable parameters. 



II. AIR SHOWER EXPERIMENTS 

Air shower detectors fall into several categories de- 
pending on the type of sensors used and on the altitude 
of the array. Scintillators such as those used in the Akeno 
array detect charge particles, which are mostly electrons 
and positrons with a fraction of order 10% of muons. In 
some arrays, the muon component can be distinguished 
from the electromagnetic component, either by a second 
layer of scintillator with an absorber between the top and 
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Ari'cLy 


Overburden (§/ cm 


type 


Energy range 


Tibet [22] 


559 


Scintillator 


1 - 200 PeV 


Akeno [15] 


909 


Scintillator/muon 


0.5 PeV - 5 EeV 


AGASA [31] 


909 


Scintillator 


4-200 EeV 


HEGRA [21] 


755 


Air Cherenkov 


0.5-10 PeV 


TUNKA [23] 


938 


Air Cherenkov 


7-1000 PeV 


Kascade [17] 


1022 


Scintillator/muon 


2-90 PeV 


Kascade- Grande [38] 


1022 


Scintillator / muon 


1-1400 PeV 


GAMMA [37] 


700 


Scintillator/muon 


3-200 PeV 


IceTop [19] 


680 


Ice Cherenkov 


1-1000 PeV 


IceCubc [54] 


680 


Ice Cherenkov (surface) + deep muons 


1-1000 PeV 


CASA-MIA [16| 


860 


Scintillator + muon counters 


0.1-10 PeV 


Fly's Eye [24] 


860 


Air Fluorescence 


1-100 EeV 


Hi-Res [29] 


845 


Air Fluorescence 


0.2-100 EeV 


Telescope Array [34] 


845 


Hybrid 


2-140 EeV 


Auger [18] 


845 


Hybrid 


1-280 EeV 



TABLE I: List of selected air shower detectors. Hybrid refers to air fluorescence telescopes looking over a surface array. 



bottom layers (as at Akeno [15J) or by separate, shielded, 
muon detectors (as in CASA-MIA [ij] and Kascade fl7j). 
The surface array at Auger [l8| and IceTop both use 
tanks filled with water or ice to detect the Cherenkov ra- 
diation produced by charged particles in the tanks. The 
method, pioneered in the Haver ah Park array [2(j, is sen- 
sitive to photons in the shower that convert in the tanks 
as well as to the less numerous e ± and muons. 



HEGRA [2l| and TUNKA [23J both use unshielded 
photomultipliers looking up at the night sky to detect 
the atmospheric Cherenkov light generated by the atmo- 
spheric cascade. In these arrays the depth of maximum, 
and hence the chemical composition of the cosmic rays, 
is judged by the lateral distribution of the Cherenkov 
light at the ground. The atmospheric fluorescence tech- 
nique, first successfully used by Fly's Eye [24] , traces the 
isotropically emitted nitrogen fluorescence lines excited 
by the passage of charged particles through the atmo- 
sphere. This method, after correcting for propagation 
through the atmosphere, maps the longitudinal devel- 
opment of each shower in the atmosphere. It therefore 
comes closest to providing a direct, calorimetric measure- 
ment of the energy of each event. By tracing each shower 
profile, a depth of maximum is assigned to each event. 
The distribution of depths of maximum is sensitive to 
primary composition. The atmospheric Cherenkov de- 
tectors provide depth of maximum and composition in- 
formation in the region of the knee up to about 100 PeV, 
while the Fly's Eye type detectors provide this informa- 
tion around 1 EeV and above. 



A partial list of air shower detectors and information 
about them is given in Table |U An account of the history 
of the air shower method and the development of the 
various techniques is given in the review of Kampert & 
Watson [25[ . The current status of spectrum and compo- 
sition is reviewed in the paper of Kampert & Unger [26| . 



III. FINDING A SINGLE ALL-PARTICLE 
SPECTRUM 



The technique of using a feature of the energy spec- 
trum to inter-calibrate different measurements was em- 
phasized by Berezinsky, Gazizov & Grigorieva in their 
model of the ankle [27J . They explain the ankle (which 
appears as a "dip" when the spectrum is plotted as 
E 3 dN/dE) as being the result of energy loss by protons 
from sources at cosmological distances to e ± pair produc- 
tion in the cosmic microwave background (CMB). They 
compared measurements of Yakutsk [28], a preliminary 
version of HiRes data [29[ , and a combination of Akeno 
(high energy) [30( and AGASA [3l|, shifting the energy 
scales respectively by 0.73, 1.2 and 0.9. In a recent re- 
view [HJ], this approach is refined, extended to include 
more recent data [33j,[34|, and compared to three differ- 
ent models of the transition from Galactic to extragalac- 
tic cosmic rays. At lower energy, the position of the knee 
has been used to provide relative calibrations of differ- 
ent measurements in the energy range from 3 x 10 14 to 
~ 10 16 eV HI. 

In this section, we show air shower data in two overlap- 
ping energy regions, 10 14 to 10 18 eV and 10 16 to 10 21 eV. 
The first includes the knee and the second the ankle 
and the end of the spectrum. In these plots the lines 
show the all particle spectrum proposed in Ref. [36| . 
which we discuss in the following section. It is impor- 
tant to keep in mind the amplification that occurs when 
data are plotted as E 1+1 dN / '&E , as is commonly done 
to display relatively small deviations from the under- 
lying power law structure of the cosmic-ray spectrum. 
Air shower experiments typically measure and present 
there results as number of events per logarithmic bin 
of energy. Uncertainties in reconstructed energy gener- 
ally scale with energy, so the resolution is expressed as 
SE/E = 5\a(E). If the energy scale is shifted by an 
amount Sln(E), then each point moves on a log- log plot 
by a distance \/T+ 7 2 S ha.(E) at an angle elevated from 
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FIG. 2: Data from kilometer-scale air shower experiments for the spectrum up to one EeV, including the knee of the spectrum. 
Left: data as presented; Right: data replotted with energies shifted as shown in the labels. 



the horizontal by 9 — tan _1 (7). For 7 = 1.6 as in Fig. [5] 
the shift is amplified by 1.89 and for 7 = 2 as in Fig. [3] 
the factor is 2.24. Thus, differences among data sets are 
not as large as they appear when multiplied by a power 
of the energy. 



A. The knee region 

To study the knee region we select several measure- 
ments that have similar (though not identical) structure 
to each other, but are offset to some extent. The left 
panel of Fig. [2] shows measurements of CASA-MIA [l|| , 
KASCADE [nl, HEGRA (nj and Tibet H3, all of which 
show a bending corresponding to the knee of the spec- 
trum. Three higher energy measurements, which start 
above the knee and extend to 10 9 GeV (Tunka [23j j . 
GAMMA [13 and KASCADE-Grande [H) are also in- 
cluded. Tibet, being a closely spaced array at high alti- 
tude has data down to 100 TeV. CASA-MIA was also a 
closely spaced array at intermediate altitude. 

The right panel of Fig. [2] shows the spectra replotted 
with energies shifted by factors as noted in the labels on 
the plot. Some differences in shape remain (particularly 
between CASA-MIA and Tibet in the knee region), but 
it is possible with moderate shifts to bring the data into 
substantially better agreement. Since the location of the 
knee is not fixed a priori, it is not possible to decide which 
experiment has the best absolute calibration. 



B. Ankle region 

Data up to the highest energies are collected in Fig. [3] 
Data in the region of the ankle are from the giant 
air shower detectors AGASA pjl, H iRes 1 and 2 [29j . 
Auger 33] and Telescope Array psij ]. The data below 
10 9 GeV from Tunka, GAMMA and KASCADE-Grande 
arc also included. The line shown is the model of Refer- 
ence [36[ with an extragalactic proton component. Larger 
shifts are needed to bring the different measurements 
into better agreement than in the case of the knee re- 
gion. In their recent review [12], the Auger data has 
been shifted up by a factor 1.22, which we also use. This 
brings the data into agreement in the decade between 10 9 
and 10 10 GeV. 



IV. APPROACHES TO DESCRIBING THE 
SPECTRUM 

Although it is possible to make a model in which the 
entire observed cosmic-ray spectrum comes from sources 
in our galaxy [3{| , it is generally accepted that the knee is 
in some way associated with the beginning of the end of a 
population particles accelerated by sources in the Milky 
Way and that the highest energy particles are from ex- 
tragalactic sources. In this section we explore two rather 
different realizations of this basic idea. In both cases we 
follow the reasoning of Peters [13] by assuming that the 
knee and other features of the primary spectrum depend 
on magnetic rigidity as defined in Eq. [T] The motiva- 
tion for this assumption is that both acceleration and 
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FIG. 3: Data from giant air shower detectors. Left: Data from AGASA 31], HiRes [29( ] . Auger pl| and Telescope Array [34| ; 
Right: same with data of AGASA shifted down in energy by 0.7 and Auger shifted up in energy by 1.22. 



propagation in models that involve collisionless diffusion 
in magnetized plasmas depend only on rigidity. The first 
evidence for a Peters cycle associated with the knee of the 
cosmic-ray spectrum comes from the unfolding analysis 
of measurements of the ratio of low-energy muons to elec- 
trons at the sea level with the KASCADE detector fl7j . 
They found that the knee occurred earlier for protons 
and helium and later for heavier nuclei. The same Pe- 
ters cycle pattern seems to occur also in the hardening of 
spectrum observed recently around 200 GV as reported 
in Refs.[| and 0. 

A. Hillas model 

The model of Ref. [36[ is an attempt to implement the 
model of Hillas ji(J in which the knee represents the end 
of the spectrum of cosmic rays accelerated by supernova 
remnants in the Milky Way and the ankle represents the 
transition to particles from extra-galactic sources. This 
picture depends on the amplification of magnetic fields by 
the turbulence associated with non-linear diffusive shock 
acceleration [4l[. Support for the presence of magnetic 
field amplification by a factor of 100 above the level the 
interstellar medium comes from the narrow rims of syn- 
chrotron radiation by electrons observed at the edges of 
some SNR [42j ]. With fields of order 100 /LtGauss, accel- 
eration of protons to energies E maK ~ 3 x 10 6 GeV is 
possible given the size and expansion rate of SNR [43| . 
In this situation it is natural to associate the knee with 
the maximum energy for the bulk of the galactic cosmic 
rays. 

If the ankle signals the transition to extragalactic cos- 





p He CNO Mg-Si Fe 


Pop. 1: 

R c = 4 PV 


7860 3550 2200 1430 2120 
1.66 1 1.58 1.63 1.67 1.63 


Pop. 2: 

R c = 30 PV 


20 20 13.4 13.4 13.4 
1.4 1.4 1.4 1.4 1.4 


Pop. 3: 

R c = 2 EV 


1.7 1.7 1.14 1.14 1.14 
1.4 1.4 1.4 1.4 1.4 


Pop. 3(*): 
R c = 60 EV 


200 0.0 0.0 0.0 0.0 
1.6 



TABLE II: Cutoffs, normalization constants (dij) and inte- 
gral spectral indexes (jij) for Eq. [3] for the implementation 
of the Hillas model (H3a) in which all populations are mixed. 
In the bottom part of the table population 3(*) consists of 
protons only (H4a). 



mic rays, and the cutoff for the SNR component occurs at 
a rigidity of several PV, then there is a gap between the 
knee and the ankle that has to be filled in by a higher en- 
ergy galactic component, which Hillas calls "component 
B." In this case there would be at least three populations 
of particles. There could of course be many more compo- 
nents in a more realistic picture in which different classes 
of sources, or even individual sources have different indi- 
vidual characteristics. For this reason a three population 
model is a minimal assumption in case the transition to 
extra-galactic cosmic rays occurs at the ankle. 

This three population picture is implemented in the 
model of Ref. [36[ by assuming that each of the three 
components (J) contains all five groups of nuclei and cuts 
off exponentially at a characteristic rigidity R c j- Thus 
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FIG. 4: Overview of the spectrum from below the knee to the ankle with the fit of Table IIIII Air shower data shifted as in 
Figs. [2] and [3] Left: lines showing individual groups of nuclei from all populations compared to data from PAMELA [s| and 
CREAM [3] at low energy. Right: shaded regions show the overlapping contributions of the three populations. 



the all-particle spectrum is given by 



ME) 



x exp 



E 



ZiR, 



c,3 



(3) 



The spectral indices for each group and the normaliza- 
tions are given explicitly in Table [TTJ The parameters for 
Population 1 are from Refs. @, Q, which we assume can 
be extrapolated to a rigidity of 4 PV to describe the knee. 
In Eq. [3]</>i is dN/dlnE and 7$ is the integral spectral in- 
dex. The subscript £ = 1,5 runs over the standard five 
groups (p, He, CNO, Mg-Si and Fe), and the all-particle 
spectrum is the sum of the five. This model is plotted as 
the solid line in Figs. [2] and [3] 



B. An alternative picture and global fit 

Spectra for the second fit are given by the same Eq. [3] 
but with qualitatively different parameters, as given in 
Table IIIII In particular, the first population has a much 
lower cutoff of R c = 120 TV. This description is related 
to the significantly harder spectra assumed for the first 
population. Each component in the first population is fit- 
ted only above R c = 200 GV, after the spectra hardening 
noted in Refs. Q and @. With these harder spectra (as 
compared to Table [TTJ) , the heavy components cannot be 
extended past the knee region. It is interesting to note 
that R c w 100TV is the classical result for the expected 
maximum energy of supernova remnants expanding into 
the interstellar medium with an un-amplified magnetic 
field of a few /xGauss [44[ . 



The spectrum with the parameters of Table IIIII is 
shown in Fig. @] from below the knee to the ankle. The 
contributions of individual groups of nuclei are shown, 
as well as the spectra of nuclei from CREAM @. We 
note that the bump in the spectrum around 10 17 eV cor- 
responds with the "iron knee" reported by KASCADE- 
Grande in their electron rich sample [45| and also noted 
by GAMMA 37|. A tendency for increasing mass above 
the knee has been noted for a long time (for example by 
CASA-MIA [46| ), which seems now to be confirmed with 
higher resolution. 

Another noteworthy feature is the possibility illus- 
trated in this fit of explaining the ankle as a Peters cy- 
cle containing only protons and iron. This possibility is 
also suggested in Ref. [HJ as an example of their "disap- 
pointing" model [47j of the end of the cosmic-ray spec- 
trum. Such a picture is disappointing because the end of 
the spectrum would correspond to the highest energy to 
which cosmic-ray acceleration is possible, rather than to 
the Greisen-Zatsepin-Kuz'min effect in which higher en- 
ergy particles lose energy in interactions with the cosmic 
microwave background |48l . l49j . 



C. Comments on fitting with several populations 

In both fits above we refer to three populations of par- 
ticles, with spectral indices for each nuclear component 
and a single characteristic maximum rigidity for each 
population. The latter assumption has the effect of mak- 
ing the composition become heavier as each population 
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p He C Fe 50 < Z < 56 78 < Z < 82 


Pop. 1: 

R c = 120 TV 


7000 3200 100 130 60 
1.66 1 1.58 1.4 1.4 1.3 


Pop. 2: 
tic = 4 rv 


150 65 6 7 2.3 0.1 0.4 

1/1 1Q1Q1Q11 1Q 1 O 

1.4 l.o l.o l.o l.Z L.Z l,z 


Pop. 3: 

rtc = 1.6 ejV 


14 0.025 

1/1 1 o 

1.4 l.Z 


Pop. 2*: 

r? a pv 

tic — 4 r V 


150 65 6 7 2.1 0.1 0.53 

1/1 10 10 
1.4 l.o l.o l.o l.Z l.Z l.Z 


Pop. 3*: 

Rc = 1.5 EV 


12 0.011 
1.4 1.2 


Pop. 4*: 
# c = 40 EV 


1.2 
1.4 



TABLE III: Global Fit results for the cutoffs, normalization constants (a^j) and integral spectral indexes (-fij) for Eq. [3] In 
the bottom part of the table(*) populations 2 and 3 are slightly modified to accommodate a Population 4 of protons to bring 
< ln(A) > down to the observed level in Fig. [5] 



approaches its maximum, as illustrated in the left panel 
of Fig. [5] Another important point is that the higher en- 
ergy populations can contribute significantly to the flux 
in the region dominated by the lower population. The 
right panel of Fig. [4] shows the overlap of the three pop- 
ulations of the global fit of Table IIIII 

The hardening of the spectrum observed by PAMELA 
and CREAM around 200 GV is suggestive of the onset of 
a new population [Ho[ • In this interpretation, the Popula- 
tion 1 of our global fit would be a higher energy popula- 
tion which becomes dominant above 200 GV, but which 
still contributes significantly at lower energies. Other ex- 
planations have been suggested. For example, Ref. 51 1 



suggests that the hardening reflects the concave spec- 
trum characteristic of non-linear diffusive shock acceler- 
ation. In Ref. [52] it is suggested that a dispersion in 
the injection spectra of different SNR is responsible for 
the hardening of the spectrum. Reference {53] shows how 
the hardening of the spectrum could be attributed to a 
change in the type of turbulence responsible for diffusion 
of the cosmic rays. 

A general feature illustrated by the various parame- 
terizations discussed here is that a Peters cycle of cutoffs 
of elemental components with rather hard spectra before 
the cutoff can produce regions of the all-particle spectrum 
that can be described approximately by steeper power 
laws. The differential spectral index between 100 GeV 
and one PeV is close to 2.6 while the index above the 
second knee, between 2 x 10 18 and 5 x 10 19 eV is approx- 
imately 3.35. The individual spectra in the global fit 
of Table IIIII for example, have differential indices below 
their cutoffs ranging from 2.2 to 2.4 (except for hydrogen 
and helium below 200 GV). 

In the case of the ankle structure, there is one model 
in which the absolute energy of the feature is fixed by 
the physical assumptions of the model. That is the orig- 
inal work of Berezinsky et al. [27| . which explains the 
dip in the plot of E 3 dN/dE as a consequence of physi- 
cal process of pair production by protons during propa- 



gation through the cosmic microwave background radia- 
tion (CMB), which fixes the energy scale. In this "dip" 
model, the extragalactic spectrum extends below the an- 
kle and the galactic-extragalactic transition occurs below 
one EeV. In this case, according to Ref. 32] there is no 
need for a second, higher energy galactic component B. 
However, in order to avoid a gap in the energy spectrum 
around 10 17 eV, the knee population would have to ex- 
tend to significantly higher the R c « 4 GV as in both fits 
in this paper. 

The different populations of particles presumably cor- 
respond to different classes of sources. For this reason 
it is instructive to compare the energy content of the 
different populations of particles and estimate the power 
required at the source. As is well known, the total energy 
in the cosmic-ray spectrum of galactic cosmic rays, which 
is dominated by particles with energy below a TeV, can 
be provided by supernova explosions at the rate of 3 per 
century. The assumption is that approximately 10% of 
the kinetic energy released goes into acceleration of cos- 
mic rays, presumably by non-linear, first order, diffusive 
shock acceleration. With 10 51 erg in kinetic energy of 
the ejecta per supernova explosion, the total power into 
cosmic rays is then ~ 3 x 10 50 erg/century or 10 41 erg/s. 

It is interesting to compare the power requirement for 
the second galactic population in the two models de- 
scribed above with the total power of the galactic cosmic- 
ray sources. To estimate this from the parameterizations 
of Tables |TT] and IIIII we start from a simplified version of 
the diffusion equation, 



N(E) = Q(E) x t bsc (E). 



(4) 



Here N(E) is the density of cosmic-ray particles (dif- 
ferential in energy) and Q{E) is the number of parti- 
cles per second per unit volume at which the sources 
inject particles of energy E. t csc (E) = t E~ s is the 
energy dependent escape time from the galaxy. We as- 
sume To = 10 7 yrs and 5 = 0.33. [62[ Multiplying by 
the factor 47r/c, which converts flux to density, we can 
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FIG. 5: Left: Mean ln(A) for the four parameterizations of tables HT1 and I lIII (For iron hi(A) ~ 4.) Right: A modified fit with 
the addition of a 4th population of extra-galactic protons (see text for discussion). 



then integrate the spectra for Population 2 given the pa- 
rameters in the Tables. We find 6 x 10 48 erg/century 
for Component B (Population 2) of the Hillas model as 
parameterized in Table HI] and 2 x 10 49 erg/century for 
Population 2 of the global fit of Table [TTTl Both are rea- 
sonable in the sense that they could be accounted for by 
special sources at a level less than 10% of the total power 
requirement for all galactic cosmic rays. 



D. Composition 

Determination of composition with indirect air shower 
experiments is difficult, primarily because of the large 
fluctuations from shower to shower, which tend to smear 
out differences arising from the different mixture of pri- 
maries. Three data samples are shown in the left panel of 
Fig.[S]below 10 18 eV. TUNKA results [H[ are from mea- 
surement of the lateral distribution of Cherenkov light on 
the ground. IceCube results are from an analysis of one 
month of data on coincident events obtained when the 
detector was half complete (54|. Coincident events are 
those in which nearly vertical showers are reconstructed 
both by IceTop on the surface and by the deep detectors 
of IceCube more than 1.5 km below the surface. The 
composition-sensitive parameter is the ratio of the ~TeV 
muons in the shower core to the shower size at the sur- 
face. Heavier primaries produce more light due to muon 
energy loss in the deep detector for a given shower size 
at the surface compared to light primaries. The KAS- 
CADE result [H is based on measurements of the ratio 
of electrons to muons in showers at the surface, and we 



take the < ln(A) > values as plotted in Ref. |26f . 

In their review (26j Kampert & Unger interpolated the 
depth of maximum measurements of several air fluores- 
cence detectors between predictions for a pure proton as- 
sumption and a pure iron assumption in order to obtain 
a value for < ln(A) >. We show in the left panel of Fig. [5] 
the values they inferred from HiRes-MIA [55| , Auger and 
Telescope array. The results depend to some extent on 
the hadronic interaction model used to calculate depth 
of maximum for protons and iron, but the trend of the 
data is similar in different models. We plot their results 
for the SIBYLL model 

An important early result was obtained by the proto- 
type of the HiRes fluorescence telescope observing show- 
ers in coincidence with the underground muon array that 
formed the muon detector for CASA-MIA. The combined 
hybrid detector allowed the profile of nearby, relatively 
small showers to be reconstructed. The measured depth 
of maximum distribution they observed [55j j showed the 
depth of maximum increasing rapidly from 10 17 eV to 
3 x 10 18 eV in a manner consistent with all the param- 
eterizations in the left panel of Fig. [S] except the disap- 
pointing model with only iron at the highest energies. 

There is a significant disagreement in interpretation of 
depth of maximum measurements above 10 18 eV, with 
Hi-Res and TA preferring nearly pure protons and Auger 
preferring a transition to heavies. But in both cases, a 
transition to pure iron, as would be the case at the end 
of the spectrum in the case of the disappointing model, 
does not seem to be indicated. The right panel of Fig. [5] 
shows a fit that contains a 4th population made entirely 
of protons to show the amount that would be needed to 
bring the composition closer to the data in this energy 
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region. There would need to comparable amounts of pro- 
tons and iron. We note that such a mixture would not be 
consistent with the small fluctuations seen by Auger , 
which would seem to require a more pure composition of 
heavy nuclei. It appears that there is not at present a 
satisfactory understanding of the highest energy cosmic 
rays. 

V. CONCLUSION 
A. Transition from galactic to extragalactic sources 

The energy range were Populations 2 and 3 meet is usu- 
ally treated as the transition from galactic to extragalac- 
tic cosmic ray sources. There are at least three different 
ways in which this transition region can be treated. 

In the dip model [27j the transition region is below 
10 18 eV. The dip itself is generated by the energy loss 
on production of electron-positron pairs in interactions 
on the microwave background. The energy spectrum of 
the extragalactic protons is flatter below 10 18 eV. What 
that means is the transition from heavy galactic nuclei to 
extragalactic protons and He nuclei happens below that 
energy. This is very different from the classical 'ankle' 
model [HI where the transition region is above that en- 
ergy and the galactic cosmic ray contribution above 10 19 
eV is significant. The chemical composition of cosmic 
rays then starts changing significantly only above 10 19 
eV. 

The third model is that of mixed cosmic ray compo- 
sition at the highest energies [Hj]. In such models the 
extragalactic cosmic rays are accelerated with a chemi- 
cal composition similar to that of the GeV galactic cos- 
mic rays. Their composition is changed in propagation 
because of photoproduction, electron positron pair cre- 
ation, and nuclear photodisintegration. The transition 
between galactic and extragalactic cosmic rays is at in- 
termediate energy, and the final energy spectrum and 
chemical composition depend strongly on the distance to 
the extragalactic sources. 

B. Outlook 

As noted in the previous paragraphs, there is much 
interest in studying the transition from galactic to ex- 
tragalactic cosmic rays, which is expected in the energy 
range between 10 17 and 10 19 eV. There is a lack of data 
at present because this energy range is at the upper 
end of the range of kilometer-scale experiments such as 
TUNKA, KASCADE-Grande and IceTop, and it over- 
laps the threshold region for the giant air shower detec- 
tors, Auger and Telescope Array. As a consequence, we 



do not yet have a good understanding of this important 
transition region. Only preliminary results from IceCube 
and IceTop have been presented so far. With the comple- 
tion of the detector and operation with the full detector 
since May, 2011, we can expect to see significant new 
results from IceCube at the beginning of the transition 
region, but the kilometer-scale detectors will always be 
statistically limited above the EeV energy range. From 
the higher energy side, the big experiments are working 
to decrease their energy threshold and become efficient 
down to the 10 17 eV energy range. 

The Auger Southern Observatory has added the High 
Elevation Auger Telescope (HEAT, jg^l) and the Auger 
Muon and Infill for the Ground Array (AMIGA, [H3) to 
the original design in order to lower their threshold. Air 
showers of energy around 10 17 eV do not emit sufficient 
light to be seen further than a few kilometers away. At 
that distance the depth of maximum X max appears at 
higher elevation, outside of the viewing angle of the orig- 
inal design. In addition, nearby showers of lower energy 
reach their maximum of development above the 30° limit 
of the original Auger telescopes. HEAT is composed of 
three fluorescence telescopes of the same basic design as 
the original Auger telescope and is installed at the west- 
ern fluorescence detector site of the observatory. They 
can operate in two positions. Horizontally they share the 
same field of view as the original telescopes. This posi- 
tion is used for calibration of the instruments. Tilted 
upward by 29°, which is the normal operation mode for 
HEAT, they observe the upper part of the atmosphere. 
The first light from one of those telescopes was seen in 
January 2009. 

An infill array of 85 detectors is deployed on two grids 
of one half (750 m) and one fourth (433 m) of the regular 
Auger surface array grid over the field of view of HEAT. 
The 750 m infill covers an area of 23.5 km 2 and the 433 
m infill covers 5.9 km 2 . There are also plans to have 
muon counters in AMIGA. The close spacing of the new 
subarray is 10 17 eV. 

The Telescope Array experiment is in the process of 
building its low energy extension TALE (G. Thomson, 
private communication). It will consist of ten fluores- 
cent telescopes located at the Northern side of the array. 
These telescopes will be elevated to be sensitive from 31° 
to 59° above the horizon. One hundred new scintillator 
counters will be deployed nearby at distances of 400 m 
from each other (1/3 of the standard spacing in the Tele- 
scope Array) . Thirty- five of them will be operational in 
March 2013. TALE will be sensitive to primary particles 
of energy above 10 16 5 eV. 

In summary, a thorough study with good statistics of 
the energy region between 10 17 and 10 19 eV should be 
scientifically productive. 
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